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I. INTRODUCTION 


The discovery of the 125 GeV Higgs boson at the LHC 7 © 8 TeV runs manifests the 
Higgs mechanism [M] and its role playing in the electroweak symmetry breaking (EWSB). 
In many of new physics models beyond the SM (BSM), the Higgs sector is extended with 
several scalar multiplets. Examples include the minimal supersymmetric extension of the 
SM (MSSM) [5], the left-right symmetric models [6], and the composite Higgs models [7j. 
There are several Higgs bosons in these models, where at least one of them is identified as 
the 125 GeV Higgs boson discovered at LHC. Therefore, those extra heavy Higgs bosons are 
yet to be searched for by the LHC experiments running at ^/s = (13 TeV, 14 TeV) and the 
future high-energy pp colliders running at ^/s = 50 — 100 TeV [H |9]. 

The two-Higgs-doublet model (2HDM) was motivated to provide extra CP-violation 
sources in the scalar sector |T0]. Such a setup is also required by the MSSM, due to the 
cancellation of gauge anomaly. To search for the second Higgs doublet, it is important to 
find the most dominant final states that are consistent with the LHC discovery and measure¬ 
ments of the 125 GeV Higgs boson. The global fits to the CP-conserving 2HDM parameter 
space were previously carried out by several groups mHHj, all of which pointed to the so- 
called “alignment limit”, i.e., Cjs-a —>• 0. In these fits, one always assumes the light CP-even 
Higgs boson h in 2HDM being the 125 GeV SM-like Higgs boson discovered by LHC. Another 
2HDM parameter to control the size of Higgs boson couplings is the ratio of two Higgs vacuum 
expectation values (VEVs) tp, with the definition given in Eq. Q. Obviously, it is impor¬ 
tant to focus on the parameter space that is allowed by the current global fits for the future 
hunting of the other heavy Higgs bosons in 2HDM. There have been considerable works on 
the LHC searches for the heavy Higgs bosons in 2HDM through different exotic decay modes, 
including A hZ/HZ [IMB], H ^ hh [HllIMS], A/H [21!, P ^ P+P’ [Mj, 

and P^ —)• AW^/HW^ [25] . The current experimental searches at the LHC 8 TeV runs 
include A —)• hZ [26ti28] and H ^ hh |26l [29l [SO] . Some of the decay modes, such as A —)■ hZ 
and H ^ hh, are due to couplings that are proportional to the alignment parameter cp-a- 
In the exact alignment limit of cp-a = 0, these decay modes will be vanishing and are of 
minor interest for the next-step experimental searches. Other decay modes, such as A ^ HZ 
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and H —?■ which involve another undiscovered heavy scalar in the final states, are 

usually suppressed by the small phase space. The conventional experimental searches for the 
heavy Higgs bosons were motivated by the MSSM scenario where the Higgs sector is dehne 
by the 2HDM-H setup. For the large-t^g inputs, the heavy Higgs boson couplings with the 
charged leptons and the down-type quarks are enhanced. Correspondingly, the important 
experimental search modes are A/H —)> {bb [31fl39j and —)• {tb^T^v) mm, 

which mostly exclude the heavy Higgs bosons from the large-t^ parameter space. On the 
other hand, the final states of ti from the heavy neutral Higgs boson decays can be quite 
important with the low- and intermediate-t^ inputs. The searches for the tt final states from 
the heavy Higgs boson decays are thought to be very challenging, where one has to deal with 
the large SM background of pp —)• tt. For the A/H production, it is known that the signal 
channel oi gg ^ A/H tt strongly interferes with the SM background |42H44] and results in 
a peak-dip structure. Therefore, one can only rely on the heavy-quark associated production 
channels oi bb+A/H and tt + A/H to study the A/H ^ tt decays. The searches for processes 
involving these hnal states at the ILC, the LHC and the future 100 TeV pp colliders were 
recently studied in Refs. [15H50] . 

In this work, we study the LHC searches for the heavy neutral Higgs bosons by tagging 
the boosted top jets from their hadronic decays. The technique of tagging boosted objects 
such as SM Higgs bosons mm, top quarks [531163] . and {W, Z) gauge bosons [M] were 
previously proposed by several different groups. Some of the previous studies of the BSM new 
state searches and the measurements of the Higgs boson properties by tagging the boosted 
jets include Refs. Generally speaking, these techniques use certain jet algorithms to 

reconstruct the entire hadronic decay of the boosted objects of {Hsm A rather than 

reconstructing the individual decay products. Specifically, there are two classes of methods 
of tagging the boosted top quarks. One algorithm is called the JHUTopTagger [^, which 
requires the summation of the transverse momenta of the decayed particles to be larger than 
1 TeV. This is very challenging when one is interested in mother particles with masses 
of several hundred GeV to 0(1) TeV. Alternatively, we study the LHC searches for the 
A/H —> tt decays by using the HEPTopTagger method, which is efficient in tagging the 
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top jets with intermediate transverse momenta of 0(100) GeV. We consider the tt + A/H 
production channels, with the sequential decay mode of A/H —)• ti. At least one hadronically 
decayed top quark will be reconstructed by the HEPTopTagger, and we denote such a top 
quark as th- For the tt + A/H channel, we look for the signal channels of th plus the same-sign 
dileptons (SSDL). 

Our paper is organized as follows. In Sec. we have a review of the heavy neutral 
Higgs bosons in the framework of the CP-conserving general 2HDM. By identifying the light 
CP-even Higgs boson h being the 125 GeV Higgs boson discovered at the LHC, the current 
global ht to the Higgs signal strengths point to the alignment limit of c^_q, —)• 0 in the 
general 2HDM. In this limit, we show that the decay modes of A/H —?■ tt are always the 
most dominant ones for 2HDM-I and also dominant ones for 2HDM-H with low-t^ inputs. 
The inclusive production cross sections of cr[pp ti+ {A/H — tt)] at the LHC 14 TeV runs 


are evaluated. In Sec. HI, we perform the Monte Carlo (MC) simulations and the top tagging 
analysis of the A/H —)• tt decay modes. We consider the heavy neutral Higgs bosons A/H in 
the mass range of G (350 GeV , 1200 GeV). One of the top quarks from the A/H —?■ tt 

decay channel will be tagged through its hadronic decay products by the HEPTopTagger 
method. For the tt + A/H associated production channel with the A/H —)• ti decays, we look 
for signal events with a boosted top jet plus the SSDL. After imposing the preselections to 
both signals and SM background, we feed the events into the ROOT TMVA package [69] to 
gain more discrimination power between the signal and SM background processes. We obtain 
a signal reach for ~ 0(1) TeV at the 14 TeV high-luminosity LHC (HL LHC) runs 

with integrated luminosities of 3 ab“^. The results are projected to the {Mj^i^j ,t^) plane of 
the general 2HDM with the alignment limit. This search turns out to be most sensitive to 
the low-ty 3 parameter space of the general 2HDM. Finally, our conclusion and discussion are 
in Sec. HVj 


II. THE HEAVY NEUTRAL HIGGS BOSONS IN 2HDM 


In this section, we briefly discuss the productions and decays of the heavy neutral Higgs 
bosons A and H in the context of the general CP-conserving 2HDM. We always take the 
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alignment limit of Cf^-a = 0 in the general 2HDM, which is consistent with the current global 
ht when the light CP-even Higgs boson h is SM-like with mass of 125 GeV. 


A. The 2HDM couplings in the alignment limit 


The most general 2HDM Higgs potential is composed of all gauge-invariant and renor- 
malizable terms by two Higgs doublets (<hi ,<1>2) G 2_|_i of the SU(2)l x U(1)y electroweak 
gauge symmetries. For the CP-conserving case, there can be two mass terms plus seven 
quartic coupling terms with real parameters. For simplicity, we consider the soft breaking of 
a discrete Z 2 symmetry, under which two Higgs doublets transform as (d>i , $ 2 ) —?■ (—‘^ 1 , ^^ 2 ). 
The corresponding Lagrangian is expressed as 




i=l ,2 


H(^>i,d> 2 ) = mn|d>ip-Fm 22 |^> 2 |^ - + P.c.)^Aild'il^-k ^A 2 |$ 2 |‘^ 

-k A3|$ip|d>2p-k A4|d>ld>2p + ^As ($|^>2)(^>ld>2) + P-c. 


( 1 ) 

( 2 ) 


Two Higgs doublets $1 and <1>2 pick up VEVs to trigger the EWSB 



and one parametrizes the ratio of the two Higgs VEVs as 


(3) 


ta = tan/3 = — . (4) 

Vl 

The perturbative bounds of the heavy Higgs boson Yukawa couplings constrain the choices 
of t/s, which should be neither as small as 0(0.1) nor as large as 0(50). In our discussions, 
we mostly focus on the parameter regions of tp ~ 0(1). The light CP-even Higgs boson h is 
taken as the only state in the 2HDM spectra with mass of 125 GeV and its couplings with SM 
fermions and gauge bosons are controlled by two parameters of (q;,/ 3). A more convenient 
choice of the 2HDM parameter set is (c^_a ,t/ 3 ). The current global fits [HHH] by using the 
LHG 7 © 8 TeV runs to the 2HDM parameters point to the alignment limit of cp-a —?• 0. 
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2HDM-I 

2HDM-H 

CU 

-IA/3 

-IA/3 

in 

-1A/3 

A 

en 

-l/A 

A 

CU 

IA/3 

i/A 

Cd 

-1A/3 

A 


-1A/3 

A 


TABLE I: The Yukawa couplings of the SM quarks and charged leptons to the heavy CP-even Higgs 
boson H and the CP-odd Higgs boson A in 2HDM-I and 2HDM-H with the alignment limit of 
C/3— a — 0- 


In the general 2HDM, SM fermions with the same quantum numbers couple to the same 
Higgs doublet, which will avoid the tree-level flavor-changing neutral currents. For 2HDM-I, 
all SM fermions couple to one Higgs doublet (conventionally chosen to be $ 2 )- This setup 
can be achieved by assigning a discrete Z 2 symmetry under which <hi —)• —d>i. For 2HDM-H, 
the up-type quarks Ui couple to one Higgs doublet (conventionally chosen to be <1>2) and 
the down-type quarks di and the charged leptons 4 couple to the other (d>i). This can also 
be achieved by assigning a discrete Z 2 symmetry under which $1 —?■ —d>i together with 
{di ,ii) —>■ {—di , —£i). Details of the 2HDM Yukawa setups were reviewed in Ref. [70]. At 
the tree level, the heavy CP-even Higgs boson H and the CP-odd Higgs boson A couple to 
the SM fermions through the Yukawa coupling terms 

= Y1 , (5) 

/ ^ 

where / is the SM fermion, m/ is the SM fermion mass, and v = -|- = {V^ = 

246 GeV. The Higgs Yukawa couplings are simplified greatly when one sets in the alignment 
limit of c^_Q = 0. The dimensionless coupling strengths of and in the alignment 
limit are listed in Table, [^for the 2HDM-I and 2HDM-H cases, which only depend on In 
both 2HDM-I and 2HDM-H cases, the dimensionless Yukawa couplings of are inversely 
proportional to Therefore, the cross sections of tt {A/H —)■ tt) are always enhanced in 
the low-t^ regions. 
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Besides the Yukawa couplings, there are also other couplings relevant to the decays of the 
heavy neutral Higgs bosons A and H. From the 2HDM kinematic terms \D^i\‘^, one has 
Higgs-gauge couplings of G{HVV), G{AhZ/AHZ), and G{AH^W^). The G{HVV) and 
G{AhZ) couplings are vanishing in the cp-a = 0 limit. From the general 2HDM potential, 
one also has the triple Higgs couplings such as G{Hhh), G{HAA), and G{HH^H~). The 
existences of these couplings lead to exotic heavy Higgs boson search strategies. 

Though we focus on the alignment limit of c^_q, = 0, it is noted that the current global 
hts to the 125 GeV GP-even Higgs boson h in 2HDM generally allow the parameter choices 
of Cjs-a ~ 0(0.1) for 2HDM-I and c^_q, ~ 0(0.01) for 2HDM-H, respectively. The relevant 
results were obtained in Refs. It was shown in the previous discussions |18L [20] 

that some of the heavy Higgs boson decay modes of H —)• hZ and H ^ hh can become 
the leading ones, especially for < 2mt. The relevant LHC searches can be performed 

by the boosted Higgs searches plus the opposite-sign-same-flavor di-leptons [18] , and via the 
66 -I -77 final states |20|. Note that the enhancements of such decay modes may rely on some 
special choices of parameters. For example, the A —)• hZ can be signihcant when masses of 
two CP-even states are nearly degenerate, and the H ^ hh mode is mostly important when 
the corresponding Higgs cubic couplings in the 2HDM potential are enhanced. In this sense, 
the full discussions of parameter regions deviating from the alignment limit rely not only on 
the values of c^_q, but also on the input parameters of the full 2HDM mass spectrum and 
the 2HDM potentials as well. Later, we present our search results in terms of the model- 
independent cross sections and map into the , tp) plane for 2HDM-I and 2HDM-H 

under the alignment limit. The results can be rescaled to other general parameter inputs 
straightforwardly. As we show below, the decay modes of A/P ti are particularly relevant 
for the low-t ,3 parameter regions of 2HDM, which was also pointed out in Ref. |71] for the 
MSSM case. 

B. The heavy quark associated productions and decays of A and H 

In the alignment limit of = 0, the production channels of the heavy neutral Higgs 
bosons A/H at the LHC include the gluon-gluon fusion and the heavy quark associated 
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processes [Tans]. The previous studies of the gluon-gluon fusion to tt final states via the spin- 
0 resonances have suggested strong interference effects with the QCD backgrounds. Therefore, 
the heavy-quark associated productions are considered as the production channels at the 
LHC. For the 2HDM-I case, all dimensionless Yukawa couplings of the SM fermions to the 
A/H are universally proportional to l/tp. One expects the scaling of the production cross 
sections (T\pp —)• hh + A/H] 1/t^ and a\pp —)• tt + A/H] l/t| with various tp inputs. For 

the 2HDM-II case, the dimensionless Yukawa couplings scale as CX/h ^ ^‘a/h ^ 

respectively. Therefore, the contributions from the top quark annihilation processes will be 
enhanced significantly with the small-t^ inputs. Note that for heavy Higgs bosons in the 
mass range of G (350 GeV, 1200 GeV), the total decay width of T can be as large 

as 0(1) — 0(10) GeV with the small-t^ inputs. Consequently, there can be potentially large 
interference effects between the signals and the QCD background. Through the evaluations 
by MadGraph5_aMC@NLO [71], it turns out that the inclusion of the interference effects 
leads to ~ 0(10%) corrections to the signal cross sections. 


A decays 

Einal states 

Alignment limit 


A -A (t+t“ ,pApr ) 

/ 

SM fermions 




A — >• (tt, hh) 

V 


A-^hZ 

- 

Exotics 

A^ HZ 

/ 



/ 

Loops 

jZ) 

/ 


TABLE If: The classification of the CP-odd Higgs boson A decay modes in the general 2HDM. A 
checkmark (dash) indicates that the decay mode is present (absent) in the c^_q, = 0 alignment limit. 


All possible decay modes of heavy neutral Higgs boson are listed in Tables. [II| and HI for A 
and H, respectively. For our discussions, we evaluate their partial decay widths in the mass 
range of £ (350 GeV, 1200 GeV) by 2HDMG [75]. Besides of setting the alignment 

limit of cp-a = 0, we also turn off all possible exotic decay modes listed in Tables. and 


HI This is reasonable when one assumes the masses of all heavy Higgs bosons are close to 









H decays 

Final states 

Alignment limit 


H -A- (r+T“ ,g^pr) 

/ 

SM fermions 



H ^ {tt, bb) 

/ 

gauge bosons 

H ^ {WW ,ZZ) 

- 


H ^ AZ 

/ 


H 

/ 

Exotics 

H —>■ hh 

- 


H ^ AA 

- 


H H+H- 

/ 

Loops 

H ^ {gg,11,iZ) 

/ 


TABLE III: The classification of the CP-even Higgs boson H decay modes in the general 2HDM. A 
checkmark (dash) indicates that the decay mode is present (absent) in the cp-a = 0 alignment limit. 



FIG. 1: Shown is cr[pp —>■ ttA/H] x BR[A/iL —>■ tt\ (for both 2HDM-I and 2HDM-II cases) with 
Ma/h G (350 GeV, 1200 GeV) at the LHC 14 TeV runs. 

each other; Ma ~ Mh ~ Mfj±. Therefore, the only relevant decay modes for heavy neutral 
Higgs bosons are A/H —?■ // and A/H —?■ gg. The loop-induced decay branching ratios 
of 'Bt:[A/H —>■ are typically smaller than 10“®, which can be neglected. For the 

2HDM-I case, the decay branching ratios of Br[A/Ff —)• tt\ are always dominant to be 0(1) 
since all dimensionless Higgs Yukawa couplings scale as ~ l/t/ 3 - For the 2HDM-H case, the 
Bt:[A/H —>■ tt\ can be suppressed to ~ 0(0.1) with the laige-tg inputs, where the partial 
decay widths of T[A/H —)> bb] and T[A/H —)> become dominant. Combining with the 
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production cross sections evaluated by MadGraph5_aMC@NLO previously, we demonstrated 
the cross sections of cr\pp —)• ti + A/H] x Bj:[A/H —)■ tt] within the mass range of G 

(350 GeV , 1200 GeV) at the LHG 14 TeV runs in Fig.[^ As it turns out, the decay branching 
ratios of Br[A/Ff —tend to unity for both 2HDM-I and 2HDM-II with the small-ty 3 inputs 
of ~ 0(1). For this reason, we combine the cross sections of a\pp —)• ti+A/H] x Br[A/H —)• ft] 
for both 2HDM-I and 2HDM-II into one plot. 

III. THE LHC SEARCHES EOR THE HEAVY NEUTRAL HIGGS BOSONS VIA 
THE ti CHANNEL 

In this section, we analyze the LHC 14 TeV searches for the heavy neutral Higgs bosons 
A and H via the ti + A/H productions, with the sequential decay modes of A/H ^ ti. We 
always tag the boosted top jets t^ by using the HEPTopTagger method. Afterwards, we 
shall look for events including a top jet th plus SSDL. The corresponding SM background 
processes should include the final states with SSDL plus multiple jets, where a jet may be 
mistagged as the boosted th- Thus, the SM backgrounds include ti [7611771, tibb, {W^Z,ZZ) 
plus jets |7S|, and {tiW^ ,ttZ) [72]. Here, we evaluate all relevant cross sections of the SM 
background processes by MadGraph5_aMC@NLO m with the NLO QCD corrections 

a{pp ti) K. 803 pb , 
a{pp —)■ tibb) « 30 pb , 
a{pp —)• W^Z) Ri 50 pb , 

a{pp —^ ZZ) « 15 pb . (6) 

To evaluate the signal significance, we should take into account the uncertainties from the 
SM background processes. For dominant SM backgrounds of ti and W^Z, the naive es¬ 
timation of the signal significance may be made by using the theoretical uncertainties of 
{6a/a)tt Ri 10% [76l [77| and {Scr/a)iY±z ~ 5% [78] due to the factorization scale uncer¬ 
tainties. Furthermore, the MG uncertainties at the extremely constrained regions of phase 
space can be typically larger. To account for the MG uncertainties, we use the uncertainties 
of ~ 20% for the ti process and the uncertainties of ~ 10% for the W^Z as the conserva- 
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tive estimations of the signal reaches. There are other SM background processes including 
(W^W^ ,tt + ,ti + Z) plus jets, with the corresponding cross sections less than 1 pb. 
As we show later, the dominant SM background processes after the preselections of plus 
SSDL are ti and W^Z plus jets. Therefore, we neglect all other SM background processes in 
our later analysis. After the reconstruction of the boosted th, we shall select the kinematic 
variables for the signal and background events and carry out the TMVA analysis to optimize 
the rate of S/^/B. The results of our analysis will be projected to the signal reaches at the 
HL LHC runs with integrated luminosities of 3 ab~^. 

A. The MC simulations and the top jet tagging 

For event generations of the signal processes, we use Universal FeynRules Output [80] 
simplified models with A or H being the only BSM particles. The relevant coupling terms 
to be implemented are the Yukawa couplings of Att/Htt. We generate events for both signal 
and SM background processes at the parton level by MadGraph5_aMC@NLO m, with 
the subsequent parton shower and hadronization performed by Pythia m- Afterwards, 
Delphes [82] is used for the fast detector simulations. In our simulations of both signal and 
background processes, we include up to two extra jets with the MLM matching in order to 
avoid the double counting. Our fast detector simulations follow the setup of the ATLAS 
detector. 


Ma/h 

350 GeV 

400 GeV 

500 GeV 

600 GeV 

700 GeV 

Rca 

2.3 

2.2 

2.1 

2.1 

2.0 

Ma/H 

800 GeV 

900 GeV 

1000 GeV 

1100 GeV 

1200 GeV 

Rca 

2.0 

1.9 

1.8 

1.8 

1.7 


TABLE IV: The choices of the jet cone sizes i?cA for different Mj^/h inputs via the tt + A/H (second 
entries) production channels. 


In what follows, we briefly describe the reconstruction of physical objects by the HEP- 
TopTagger method. The details of the algorithm can be found in the original Refs. |54H56j . 


II 














FIG. 2: The tagging rates of the top jet versus the CA jet cone sizes Rqa G (1.0,3.0). Two 
samples of Ma/h = 500 GeV (blue curve) and Ma/h = 1000 GeV (green curve) are shown. 


The energy flow observables from the Delphes output are used for the jet substructure anal¬ 
ysis by Fastjet [83]. In each event, we cluster the top jets by using the Cambridge-Aachen 
(CA) algorithm jHH |85] with certain jet cone size Rqa- By setting the reconstructed top 
mass range of G (140 GeV , 210 GeV), the HEPTopTagger algorithm finds a candidate 
boosted top jet which contains three subjets with their total transverse momenta greater 
than 200 GeV. The rate of tagging one th can be ~ 30% — 60% with certain choices of Rqa- 
As an illustration, we plot the tagging rates of the boosted top jet th versus the GA jet cone 
sizes for the Ma/h = 500 GeV and Ma/h = 1 TeV samples in Fig. It is also likely to 
tag two boosted top jets at the rate of ~ 10% — 20%. For such cases, we always choose the 
one with the largest pT as the th- Generally speaking, the tagging rates of top jets vary with 
different choices of the jet cone sizes Rqa- The boost factors of top jets are enhanced with 
the heavier resonances of Ma/h- For each mass of signal process pp ^ tt + A/H, we scan 
the jet cone sizes Rqa G (1-0,3.0) at the step of 0.1 for reconstructing the top jet th in the 
HEPTopTagger. In addition, the effects due to the underlying events can be eliminated by 


the filtering procedure ISD in the HEPTopTagger. In Table. |IV[ we list our choices of Rqa 
for the sample models with masses of Ma/h G (350 GeV, 1200 GeV). The selection criteria 
is given below for two different production channels. Once a th is obtained in an event, we 
erase its constituents from the input particles. The remaining particles will be clustered into 
narrow jets by using the anti-fcf algorithm with a jet cone size of i?narrow = 0.4. The narrow 
jets are required to satisfy pT > 20 GeV and \r]\ < 4.5. 
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Ma/h 

Signal 

tt 

ttbb 

ZZ 

s/Vb 

Total cross section (fb) 

50 

8.0 X 10® 

3 X 10^ 

5.0 X 10^ 1.5 X 10^ 


Preselection of th (fb) 

28 

3.1 X 10® 

1.4 X 10^ 

7.7 X 10® 1.9 X 10® 


Preselection of t/j-l-SSDL (fb) 

0.48 

0.56 

0.11 

3.92 0.17 

2.3(1.2) 


TABLE V: The preselection efficiencies of the Ma = 500 GeV (with Rqa = 2.1) signal and back¬ 
ground processes at the 14 TeV HL LHC. We assume the nominal cross section for the signal process 
to be a[pp ti+ A/H] x BR[A/H —>■ tt] = 50 fb. The signal significances are obtained by the naive 
estimation (the conservative estimation) of the SM background uncertainties. 


B. Multivariable analysis 

For the tt {A/H —>• tF) signal channel, one has four top quarks in the final states. As 
stated previously, at least one top quark th from the A/H decay should be reconstructed 
by the HEPTopTagger method through its hadronic decay mode. Afterwards, we select 
events containing SSDL from the semi-leptonic decays of two other top quarks. It turns 
out a significant suppression to the SM background can be achieved by selecting the events 
containing th plus SSDL. An example of the preselection efficiencies of events for the Ma/h = 
500 GeV case is tabulated in Table. |V} The suppression rates of SM background events from 
the tt and ttbb can be as significant as ~ 10“® when imposing the SSDL selection criterion. 
Obviously, the W^Z background becomes the most dominant one after the preselections. 
Meanwhile, one has cr(tt)seiect ~ 0.1 (T(IT^Z)geiect after the preselections. In Fig. we 
display the cross sections of W^Z and tt after selecting the th plus the SSDL events versus 
different choices of the jet cone sizes Rca- 

To further discriminate the signal events from the SM background, we adopt the MLP 
neural network analysis in the ROOT TMVA package. The list of kinematic variables for the 
multivariate analysis include {pt,V^4’) of i'f^h ,Ff ,if), 1 /t, number of (6-jets, non-6jets), 
PT{bo,jo), YhjPTii)-, and J2bPT{b)- Here jo and 6o denote the leading non-6-jet and the 
leading 6-jet ordered in their transverse momenta, respectively. A 6-tagging rate of 70 % |86] 
is assumed to reconstruct the 6-jets. For signal events with different masses of Ma/Hi Ha® Pt 
distributions of the tagged th are shown in the upper left panel of Fig. For the samples 
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FIG. 3: The SM background cross sections after the th plus SSDL preselections versus the CA 
algorithm jet cone sizes. 

shown here, the pT of the tagged th is not a very significant discriminator among different 
signal samples. Dne to this reason, we find the signal reaches for different inputs 

are pretty close within the mass range of the study. On the upper right panel of Fig. we 
display the distributions of the pr(jo) between the signal processes versus the SM backgrounds 
for the = 500 GeV signal sample. There is clear discrimination between the signal 

and background events on the pT of narrow jets. Afterwards, we feed the events passing the 
preselections into the TMVA package for further optimization of the cuts among all kinematic 
variables listed above. As it turns out, the kinematic variables oi pT{bo ,jo) and the number 
of 6-jots are the leading important ones for the TMVA to apply cuts. The discriminations 
between signal and SM backgrounds from the TMVA analysis are also presented for the 
Ma/h = 500 GeV and Ma/h = 1000 GeV samples. After obtaining the cut efficiencies, we 
convert the results to the signal cross sections within the 5cr discovery limits, which read 
a[tt + {A/H —>■ tt)] ~ 5 — 8 fb (naive estimation) or a[tt + {A/H —>■ tF)] ~ 10 — 15 fb 
(conservative estimation) with Ma/h £ (350,1200) GeV. The results are demonstrated in 
the left panel of Fig. By looking for the th plus the SSDL signals, our analysis shows that 
the HL LHC searches are likely to reach the heavy neutral Higgs boson masses up to 0(1) TeV 
in the low-t^ regions for the general CP-conserving 2HDM. The model-independent signal 
cross sections for the 5 a reaches are further projected to the {Ma/h ) plane, as shown in 
the right panel of Fig. The previous Ref. [38] studied the inclusive searches for the four 
top final states from the 2HDM heavy Higgs bosons at the LHC, where they obtained a mass 
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FIG. 4: Upper: The pr of the th for the signal processes tt + {A/H —>• tt) with Ma/h = 

(500 GeV, 800 GeV, 1000 GeV) (left panel) and the pr of the leading narrow jet reconstructed 
by the anti-fcj algorithm for Ma/h = 500 GeV (right panel) and SM background processes after the 
preselections. Lower: The normalized distributions of MLP neural network response for signal and 
background for the tt + {A/H —)• tt) channel, left: Ma/h = 500 GeV, right: Ma/h = 1000 GeV. 

reaches of 0(800) GeV with tp < 1.5 at the 2 a C.L.. As a comparison, our search strategies 
are likely to improve the signal reaches for four top final states with mass up to 0(1) TeV 
at the 5 cr C.L. thanks to the reconstruction of boosted top jets and the TMVA analysis. 


IV. CONCLUSION AND DISCUSSION 

In this paper, we have carried out an analysis of the heavy neutral Higgs boson searches 
via the tt decay modes at the LHC 14 TeV runs. In the simplified scenario where one sets 
the alignment limit of Cj^-a = 0 and turns off all possible exotic decay modes of heavy neutral 
Higgs bosons, the decay branching ratios of tt final states can be usually approaching 0(1) 
with low-t^ inputs. Correspondingly, the searches for the A/H —)• tt are of the top priority 
from the perspective of the production cross sections. We consider the tt + {A/H tt) signal 
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Ma/h (GeV) (GeV) 


FIG. 5: The signal predictions together with the signal reaches (dashed lines) of ti + {A/H —)■ tt) at 
the HL LHC runs. Upper: The mode-independent cross section reaches at the HL LHC runs. Lower: 
The 5 a (in blue) and 3 cr (in green) signal reaches projected on the (M^/^ , t^) plane at the HL LHC 
runs. We present results by using both naive estimations (left panels) and conservative estimations 
(right panels) of the SM background uncertainties. 

channels in this work, whose interference effects with the QCD background are less severe 
compared to the gluon fusion channel. In order to suppress the corresponding SM background 
contributions, we adopt the HEPTopTagger method to tag the boosted top jets for both 
signal processes. Because the boost factor thus the spreading angle of th is mainly controlled 
by the mass of A/H, for each input, the cone size parameter of the HEPTopTagger 

is optimized to gain the maximal signal top tagging rate. Eor the tt + {A/H —)■ ft) signal 
channel, we require a boosted t^ plus the SSDL in the final state. Afterwards, the MLP neural 
network analysis is applied based on the TMVA package. Eor G (350 GeV , 1200 GeV), 

the production cross sections oiti+{A/H tF) as small as ~ [5 —8] fb (naive estimation) or 
~ [10 — 15] fb (conservative estimation) can be discovered at 5cr C.L.. Based on these results, 
we eventually obtain the LHG signal reaches for the heavy neutral Higgs bosons via the tt 
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final states on the plane. Our results are most sensitive to the heavy neutral 

Higgs boson searches with the low-tp inputs, which are complementary to the previous LHC 
searches performed for the {bb final states. For the 2HDM inputs of tjs ~ 0{1) with 

the alignment limit, the < 1 TeV region can be reached at the HL LHC runs via the 

boosted top jet tagging plus the SSDL final state searches. 
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